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Microstructural characterization of synthetic periodic multilayers by x-ray standing waves have 
been presented. It has been shown that the analysis of multilayers by combined x-ray refiectometry 
(XRR) and x-ray standing wave (XSW) techniques can overcome the deficiencies of the individual 
techniques in microstructural analysis. While interface roughnesses are more accurately determined 
by the XRR technique, layer composition is more accurately determined by the XSW technique 
where an element is directly identified by its characteristic emission. These aspects have been 
explained with an example of a 20 period Pt/C multilayer. The composition of the C-layers due 
to Pt dissolution in the C-layers, Pt^Ci-^,, has been determined by the XSW technique. In the 
XSW analysis when the whole amount of Pt present in the C-layers is assumed to be within the 
broadened interface, it leads to larger interface roughness values, inconsistent with those determined 
by the XRR technique. Constraining the interface roughness values to those determined by the XRR 
technique, requires an additional amount of dissolved Pt in the C-layers to explain the Pt fluorescence 
yield excited by the standing wave field. This analysis provides the average composition Pt^Ci-x 
of the C-layers. 

PACS numbers: 07.85.-m, 61.10.Kw, 68.35Dv 



I. INTRODUCTION 

Improvements in the thin film deposition techniques 
in recent years have led to the fabrication of layered syn- 
thetic microstructures (LSM) consistingof thin layers of 
alternating elements or compounds |l||J. These mate- 
rials have unique structural ||, magnetic Q and elec- 
tronic [^| properties with a wide range of applications. 
LSM containing alternating layers of high atomic num- 
ber elements (eg., W, Mo, Pt etc.) and low atomic num- 
ber elements (eg., C, Si etc.) are being used as x-ray 
reflectors ||. Indeed, x-ray multilayer optics are now 
used in many applications including x-ray astronomy, mi- 
croscopy, spectroscopy, as filters and monochromators for 
intense sources such as synchrotron radiation and x-ray 
laser cavities. It is important to correlate the measured 
properties with structure so that preparation techniques 
can be optimized to yield high performance materials. 
X-ray techniques are very useful for the measurement 
of the microstructural aspects of the multilayered sys- 
tems. Here we present the application of combined x-ray 
standing wave and x-ray refiectometry techniques for the 
microstructural analysis of periodic multilayers. 

For a perfect single crystal, according to the dynami- 
cal theory of x-ray diffraction |]],|| , a standing wave field 
is generated within the crystal as a result of superposi- 
tion of the incident and the diffracted waves when x-rays 
are Bragg-reflected by the crystal. The equi-intcnsity 
planes of the standing wave field are parallel to and have 
the periodicity of the diffracting planes. At an angle of 
incidence corresponding to the rising edge of the diffrac- 
tion peak the antinodal planes of the standing wave field 
lie between the diffracting planes. As the angle of inci- 
dence increases the antinodal planes move continuously 
inward onto the diffracting planes at the falling edge of 
the diffraction peak. Over the angular region of Bragg 



reflection, emission such as fluorescent x-rays [p|~pd| and 
electrons Jl2| from the crystal is strongly modulated, be- 
ing maximum (minimum) when the antinodal (nodal) 
planes coincide with positions of the atoms in the crystal 
or on the surface. By measuring the angular dependence 
of the intensity of the emitted fluorescence and compar- 
ing with the computed angular dependence, the standing 
wave field has been used as a structural probe to de- 
termine the positions of the impurity atoms in crystals 
-0,[l3|K adsorbed atoms on surfaces O), atoms at a 



layer/substrate interface 15 and to study thermal effects 
such as broadening of atomic position due to thermal vi- 
bration |l6| and order-disorder transitions JT?]]. Various 
applications of the x-ray standing wave (XSW) technique 
to problems relating to single crystal surfaces and inter- 
faces may be found in recent reviews |l8|JT^ ]. 

The standing wave phenomenon was also observed in 



multilayer mirrors 
tilayer films §1. T 



23] and Langmuir Blodgett mul- 
'his standing wave field was also used 
in different ways for analyzing the local structure of mul- 
tilayers p4| , p5|l , density evaluation of deposited films on 
multilayers ]26[ | and selective extended x-ray absorption 
fine structure analysis 

For a periodic multilayer system x-ray reflectivity 
(XRR) is used to determine bilayer periodicity, interface 
roughness and the fractional thickness of the layers in 
a bilayer. Interface roughness characterization by x-ray 
standing waves has been attempted for a Ni/C multi- 
layer system j25|. However, the extracted parameters 
were not optimized. Matsusita et al. (28) used XSW 
to determine the density of impurity atoms in a multi- 
layer structure. Here we present a combined reflectivity 
and standing wave characterization of a periodic mul- 
tilayer system to extract various structural parameters. 
As an example we use a 20-pcriod Pt/C multilayer sys- 
tem. Comparing with experimental data we show that, 
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structural parameters extracted from x-ray reflectivity 
analysis cannot explain the Pt fluorescence yield excited 
by x-ray standing waves. Explanation of the Pt fluo- 
rescence yield, additionally requires the presence of an 
amount of dissolved Pt in the C-layers. XSW analysis 
provides the amount of dissolved Pt in C and the av- 
erage composition Pt^Ci-x of the C-layers. Probing a 
small quantity of material dissolved from one layer into 
the other layer of a layer-pair in a multilayer system is 
very important for magnetic multilayers where alternat- 
ing layers are magnetic and nonmagnetic materials. A 
small amount (even a few percent) of magnetic impurity 
(either from the magnetic layer or external) in the non- 
magnetic layer can change mag netic coupling and mag- 
netoresistance significantly [p29j . presumably because of 
changes in the topology of the Fermi surface of the re- 
sulting alloys.. The importance of the combined XSW 
and XRR analysis is elucidated. 



where 



(4) 



In Eqn.(4) No is Avogadro's number, p m is the mass den- 
sity of the element in the layer with atomic number Z and 
atomic weight M. /' and /" are the real (dispersive) and 
the imaginary (absorption) anomalous dispersion factors, 
respectively, p is the electron density (including disper- 
sion) and /i is the linear absorption coefficient for the 
incident photons in the medium. r e is the classical elec- 
tron radius. We consider the medium for the incident 
beam to be vacuum with eo = 1. 

For the s-polarization of the electric field and smooth 
interfaces the complex coefficient of reflection Tj and 
transmission tj, being the ratio of electric fields at the 
j, j+l interface, are given by Fresnel's formulas 



II. THEORY 

We give a brief theoretical background for the x-ray 
standing wave generation inside a multilayer system. We 
mainly follow the treatment given by Dev et al. |30| for 
the formation of standing waves and resonance enhance- 
ment of x-rays in layered materials using the recursion 
method of Parratt |51[ . Then we obtain the field inten- 
sity for a periodic multilayer system and compute the 
angular variation of fluorescence yield from constituent 
elements in the multilayers. The fluorescence yield pro- 
file depends on the structural parameters of the multi- 
layer. A consistent set of microstructural parameters of 
the multilayer is obtained from the combined analysis of 
reflectivity and fluorescence yield. 



A. Reflection from a multilayer system 

If all interfaces are parallel in a multilayer system Fig. 
[j], a plane electromagnetic wave of frequency u> in a 
medium j at a position r can be written as 



tj = 



(5) 



(6) 



For small 5j,f3j approximation, no distinction need be 
made between s-polarization and p-polarization |3f| . 

For rough surfaces these expressions are to be modi- 
fied. There are several methods for o btaini ng modified 
expressions. In a well-known method 
plied by a factor Sj given by 



53| rj is multi- 
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2 (Tj kj ; 2, k 



(7) 



where <Tj is the root-mean-square deviation of the in- 
terface atoms from the perfectly smooth condition. An 
expression like Eqn.(7) is only valid for small roughnesses 
< !)• 

For the modification of tj , it is to be multiplied by 



Tj = exp[<jj (kj 



72] 



(8) 



Ej(r) = Ej(0)exp[i(ujt - k,-.r)] 



(1) 



where Ej (0) is the field amplitude at the top of the j-th 
layer. 

For all j, the components of the wave vector, kj = 
k' — ik" , are given by 
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kj.x = —cosO ; kj tZ — —r-(tj — cos 2 9) 1 ^ 2 
A A 



(2) 



So far we have discussed the reflection and refraction 
at a single interface. For a multilayer system, involving 
multiple interfaces, the electric fields at all the interfaces 
can be obtained from either a recursion relation or from 
a matrix formalism. In the following we will use the 
method of recursion relation. In the recursion method 
|3l],|4|], the transmitted field Ej and the reflected field 
Ej at the top of the j-th layer are found from the rela- 
tions : 



where 9 is the glancing angle of incidence, A is the wave- 
length of the incident x-rays and the dielectric function 
tj is given by 



1 25j - i2f3j 



(3) 
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(9) 



(10) 



2 



and 



where 



X; 



( r 3 S j +a 2 3+i x j+i) 
1 + a j+1 X j+l rjSj 



exp(- 



■ikj >z dj) 



(11) 



(12) 



dj being the thickness of the j-th layer. For the substrate 
£[ = JQ = . 

The electric field amplitudes E* (transmitted) and EJ 
(reflected) can be computed from the knowledge of A, 9, 
ej's, the thickness of the layers (dj's) and the interface 
roughness (ct/s) using Eqn.(2) through Eqn.(12) and the 
reflectivity R is then obtained from the ratio of E-fields 
outside the surface : 



R(9) = 1*5/13*1 



(13) 



For reflectivity from a periodic synthetic multilayer sys- 
tem involving interface roughness, this treatment is es- 
sentially equivalent to that of Underwood and Barbee 

For a periodic multilayer system, below the critical 
angle of incidence, 91=^26\, there exists an evanescent 
wave below the surface and total external reflection of 
the incident beam occurs {\E r D \ w \E t Q \). The interfer- 
ence between E r Q and E a can form standing waves above 
the surface |p2| . For 9 > 9\, the incident beam pene- 
trates into the first layer of the multilayer system. When 
6 1 > $2; the incident beam penetrates into the multilayer 
system for 9 > 9\. If 9% > 9\, there is the possibility 
of resonance enhancement of x-rays in medium T' for 
91 < 9 <9\ @||,|6). For 9 greater than both 9\ and 
#2, the x-ray beam penetrates into the multilayer and if 
the multilayer is periodic, Bragg diffractions can occur 

For a periodic multilayer system of x-ray reflectors the 
multilayer period consists of one low and one high elec- 
tron density alternating layers (say, Pt/C/Pt/C ...), the 
higher density layer works as a marker and the low den- 
sity layer works as a spacer. This arrangement makes 
the system an artificial periodic structure. Therefore, in 
the reflectivity from such a periodic multilayer system, 
Bragg peaks appear at positions determined by Bragg's 
law (including refraction and absorption). 



or 



2(d 1 k[ z + d 2 l4 z ) = 2nix 



2{d\sin9i + d2Sin92) — n\ 



(14) 



(15) 



where the period of the multilayer is d = di + d2 and n 
is the order of reflection. 

It is well known from the dynamical theory of x-ray 
diffraction from perfect crystals that Q a standing wave 
field is set up in the crystal during diffraction. The antin- 
ode position of this wave changes over half the unit-cell 
distance in passing the diffraction peak. This is also true 
for x-ray diffraction from a periodic multilayer system 
which will be illustrated later. 



B. Field Intensity 



The interference between the incident E-field (Ej) and 
the reflected E-field (Ej) can form standing waves within 
any layer. In order to obtain this standing wave field in 
the j-th layer one needs to know the fields E* and Ej as 
a function of depth (z). The total E-field at a point r in 
the j-th layer is given by 



Ej{r) = EHv) + EUv) 



where 



EUy) = E t ,(Q)exp{-ikj^z)exp[i(ijjt - 



(16) 



(17) 



and 



R r (r) = EUO)exp(+ik jtZ z)exp[i(ojt - k jtX x)] (18) 



Here the origin has been chosen to be on the interface 
at the top of the j-th layer. Thus ^'(O) and -BJ(O) repre- 
sent the transmitted and the reflected E-fields at the top 
of the j-th layer. E*(0) and £J(0) are readily obtained 
from the recursion relations (Eqn.(9) through Eqn.(12)). 



The field intensity 1(9, z) 



|i^ T (r)| 2 is given 



by @ 
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1(9, z) = |4(0)| 2 [ e: rp{-2fc;>} + \-J F -±\\ X p{2kl z z} 



^J-\cos{p(9)+2k' :j;Z z}}, 



where v(9) is defined by -gr^ 
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(19) 



.e., v(9) is 



the phase of the E-field ratio at the top of the j-th layer. 
If the absorption in the medium is ignored ( i.e, k" z = 
), Eqn (19) reduces to 



/(M = |£j(0)| 2 [l + 

EUO) 
+n^\cosM9) 



| ^(°) |2 

1 W 

2kj z z}], 



(20) 



It is clear from Eqn. (19) and (20) that a standing wave 
is generated within the j-th layer. The quantity within 
the square bracket in Eqn. (20) may attain a maximum 
value of 4, for \Ej(0)/Ej(0)\ =1. For small angles of inci- 
dence (9) , in some situations there are possibilities of res- 
onance enhancement of x-ray intensity in the layer. This 
has been described in details by Dev et al. However, 
at 9 > 91 and 9%, ^ (Q) / E](Q)\ 2 < 1 for a nonperiodic 
multilayer, and the field intensity is essentially given by 
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the first term in Eq. (19) or (20) with a slight modu- 
lation from the second and the third terms. For such 9 
values the reflectivity is only significant when 8 satisfies 
the Bragg condition for reflection from a periodic multi- 
layer. Standing waves are set up in the multilayer when 
Bragg diffraction occurs. This can be seen from Eqn.(20) 
by inserting the Bragg condition [Eq. (14)] 

2{k' lz di + k' 2z d 2 ) = 2k' z d = 2nir 

or 



where k z is the weighted average value for a layer-pair 
of the multilayer with periodicity d = d\ + d 2 . While 
the magnitude of the E-field ratio varies to some extent 
for layer 1 and layer 2 of the bilayer, we can approxi- 
mate this to be equal to its value just above the surface, 

[i- e -» l§Sjl 2 ~ l|fl 2 = from Eqn.(13)]. Now for 

normalized incident intensity, inserting the value of k' z in 
Eqn.(20) we obtain (for n—1) 

1(9, z) = l + R(9) + 2^R(9) cos{v(9) + ^z} (22) 

It is clear that Eqn.(22) now defines a standing wave 
within the multilayer with a periodicity d and has the 
same form as that derived from the dynamical theory 
of x-ray diffraction from perfect crystals P,|l8|]. In the 
dynamical theory of x-ray diffraction the E-ficld in a 
medium is calculated by solving Maxwell's equations in 
that medium and obtaining solutions consistent with the 
Bragg's law. This E-field, then, describes the x-ray 
standing wave intensity as a function of angle over the 
region of the Bragg peak where the phase of %t(9), v {Q)i 
changes by 7r radian |ffl ]llJ37| . The actual value of u(9) 
on the higher-angle side beyond the diffraction peak de- 
termines the position of the diffraction planes p7j ]. In 
order to show the similarity between the expressions for 
the standing wave intensity in the dynamical theory for 
perfect crystals and in the present case for multilayers 
we have inserted the Bragg's law into Eq. (20) and ob- 
tained Eq. (22), which is the well-known form obtained 
from the dynamical theory, where T/cf is the magnitude 
of the reciprocal lattice vector for the concerned diffrac- 
tion. The phase variation, v(9), for the present case of 
multilayer is shown in Fig. |^. This has the similar form 
to that obtained from the dynamical theory Jl^j37|] . 

A periodic multilayer structure can be characterized 
by generating standing waves within the multilayer and 
measuring the standing-wave-excited fluorescence yield 
from one or more elements present in the multilayer. This 
is explained in the following sections. For the computa- 
tion of standing wave field intensity, 1(9, z), we will use 
the more rigorous form of Eq.(19). 



C. Examples of calculation 

In this section we present the results of calculations of 
various quantities in sections II. A and II. B using an ex- 
ample — a periodic multilayer system consisting of 20 bi- 
layers of Pt/C on a glass substrate. The discussions pre- 
sented here are general and not restricted to only Pt/C 
multilayers. 

For multilayers, earlier analyses were performed as- 
suming the same roughness for both types of interfaces 
(A/B and B/A) in the multilayer (A/B/A/B...) @Q. 
In general, these values should be different. Surface free 
energy of the materials, a a and <tb, partly control the 
interface morphology during the growth. If a a < <Jb, 
it is the wetting condition for the growth of material A 
on material B and a nonwetting condition for the growth 
of material B on material A. Thus A-on-B (A/B) inter- 
face is expected to be smoother. The situation would be 
reverse for a a > <Jb- Indeed, high resolution electron mi- 
croscopy on W/C multilayers shows that the interface of 
C growing on W is much sharper than that of W growing 
on C p9[ . It must be noted that aw > &c- However, 
other factors such as growth temperature and interdif- 
fusion or chemical reaction between species across the 
interface also affect the interface roughness pOf . In any 
case, there is no reason to assume the interface roughness 
for both types of interfaces to be equal. Here we assume 
different roughnesses for the Pt-on-C (a\ ) and the C-on- 
Pt (172) interfaces. It will be shown later that we indeed 
get a better fit to experimental data when o\ and o 2 are 
allowed to be different. 

In Fig. [2|, we show the simulated reflectivity curve for 
smooth surface and interfaces alongwith those for sev- 
eral sets of values of surface and interface roughness. To- 
tal external reflection at low angles and multilayer Bragg 
peaks upto fourth order are seen. The higher order peaks 
are more drastically affected by the surface (oq) and in- 
terface roughness (01,02). The spacing between Bragg 
peaks is determined by the periodicity or the bilayer 
thickness (d). So these parameters can be determined 
from the reflectivity data by a least-squares fitting pro- 
cedure. In these computations we have used ept= 1 — 
(2.302xl0~ 5 ) - i(2.596x 10~ 6 ) and e c = 1 - (3.016x 
10~ 6 ) - i(8.138x 10- 10 ), (p Pt = 5.05 electrons / A 3 , 
pc = 0.698 electrons/ A 3 ), A = 0.709 A (MoK ai x-rays) 
and d = 43 A (di = 17 A,d 2 = 26 A). X-ray standing 
wave intensities are shown in Fig. [5] over the first Bragg 
peak region (9 = 0.3° to 9 — 0.6°) at several angles shown 
on the reflectivity curve in the inset. The variation of 
phase, v(9), of -§r|jjy and -§r|{jy are shown in the second 

inset of Fig. ||. The field intensity, I(z), can be ob- 
tained using R(9) and v(9) from the insets and Eqn.(22). 
However, we have used the more rigorous Eqn.(19) to 
compute the field intensity I(z) at several values of 9. 
At an angle away from the strong reflection region (a) 
the field intensity, I(z), has a weak modulation around 
a value of unity. At the low-angle side of the diffraction 
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peak (b), there are antinodes of the standing wave field 
in the C-layers (nodes in the Pt layers). As 9 increases 
(b^c^d^e), the antinodes shift inward and finally co- 
incide with the Pt-layers. The field intensity over the 
Pt-layers gradually increases as 9 increases. The inte- 
grated field intensity in the Pt-layers, 



IptiO) 



E 

j—odd 



Ij(6, z)dz 



(23) 



is shown in Fig. ^. Ip t (9) for smooth surfaces and inter- 
faces (do = <J\ = f2 = 0) and for several sets of <to, <J\, o~2 
values also shown. It is clearly seen that the field in- 
tensity 1(9) variation with 9 is sensitive to surface and 
interface roughness. The integrated field intensity over 
the carbon layers, 



Ic(6) 



j—even 



Ij(Q, z)dz 



(24) 



for (To = ci = a 2 = is also shown in Fig. ^. It is no- 
ticed that the field intensity in the Pt-layers peaks at the 
high-angle edge while the intensity in the C-layers peaks 
at the low-angle edge of the reflectivity peak. This op- 
posite trend holds the clue for the determination of the 
concentration of any dissolved Pt in C-layers. 

Our objective is to find the Pt distribution in the Pt/C 
multilayer. In the dipole approximation, fluorescence 
yield from an atom is proportional to the field intensity 
on the atom. Thus with the measurement of fluorescence 
yield from Pt, it is possible to determine the Pt distribu- 
tion. Fluorescence yield from Pt in the Pt layers should 
follow curve T' in Fig. ^, while fluorescence yield from 
Pt in the C-layers should follow curve '2'. So the ef- 
fective shape of the fluorescence yield curve will depend 
on relative concentrations of Pt in the Pt-layers and the 
C-layers. 

Interface roughness can be due to actual roughness 
or diffusion across the interface. The Pt distribution, 
f(z), with interface roughnesses <ti 7^ 01 is schematically 
shown in Fig. [|. It is obvious that a fraction of Pt is 
in the C-layers near the interface. The fluorescence yield 
of Pt generated (I-^ 9 ) from any depth is proportional to 
the product of the field intensity and Pt concentration at 
that depth. 



if 9 1 



,z) = ci j {e,z)f j {z) 



(25) 



where C is a constant. The fluorescence yield detected 
outside the sample is given by 



f{6,z) = CI J (9,z)f J (z) x exp[-^( V d m + z)] 



m=Q 



(26) 



with do=0 and the depth integrated detected fluorescence 
yield is 



N j-i 

lV(9)=Cj2expl-^(J2d m )} 
£ — ' sma. — ' 

j — l m—0 



Ij (8, z)fj (z)exp(— ^° ut z)dz 



sma 



(27) 



where a is the angle between the sample surface and the 
direction of the fluorescence detector from the center of 
the sample surface, and (i ou t is the weighted average lin- 
ear absorption coefficient for the outgoing (fluorescent) 
photons. 

The distribution of Pt concentration over the bilayers 
across the Pt-on-C interface is given by 



(28) 



for —d\ < z < d%, z—Q is on the Pt-on-C interface, o\ is 
the Pt-on-C interface roughness. Pt distribution across 
the C-on-Pt interface is given by 



(29) 



for — di < z < d\ where 2=0 is taken on the C-on-Pt in- 
terface. cr% is the C-on-Pt interface roughness. The total 
Pt distribution f(z) over the bilayer and two interfaces 
is schematically shown in Fig. ||. f(z) — fi(z) + f%{z) 
in the C-layers whereas in the Pt layer f(z) — fi(z) or 
/2(z), whichever is lower. The interface roughnesses o\ 
and iT2 are those used in the analysis of reflectivity. Now 
that the Pt distribution, f(z), over the total thickness 
of multilayer is defined, the integrated detected fluores- 
cence yield, I^ d (9), can be computed using Eqn. (27). 
The Pt fluorescence yield computed for this distribution 
of Pt over the first order Bragg reflection angular region 
is shown in Fig. |^. 

The solid curve (<tq = 3, o\ = 5, 02 = 3 A) in Fig. 
^ shows the computed fluorescence yield profile for Pt 
only in the Pt-layers. In this calculation the effect of 
roughness enters only in the computation of field inten- 
sity and the contribution to fluorescence yield from Pt 
in the C layer due to interface broadening is neglected. 
This means, in the Eqn. (27), only the sum over j=odd 
layers have been considered. Sum over all layers contain 
the fluorescence yield contribution from Pt distributed in 
the C-layers as well. The fluorescence yield curve includ- 
ing this contribution is shown by the dashed line (/ c =l, 
the significance of f c will be discussed later). 

The possibility of a small amount of dissolved Pt in 
the C-layers, in addition to the Pt in the interface pro- 
file, has not yet been taken into account. In the com- 
putation of reflectivity the existence of such dissolved 
Pt in C should enter as a change in electron density of 
the C-layers. However, due to the low electron density 
of C (0.698 electrons/ A 3 ), reflectivity is not very sensi- 
tive even to a relatively large change in C-layers electron 
density. Reflectivity for a 15% higher electron density 
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(0.803 electrons / A 3 ) of the C-layers, shown in Fig.f7|, is 
hardly distinguishable from that for pure C electron den- 
sity. Moreover, the electron density of the C-layers not 
only depends on the amount of dissolved Pt, but also on 
the change in C-layers thickness upon Pt incorporation. 
The electron density can also change due to incorporation 
of ambient atoms (e.g. Ar) during multilayer deposition 
p8| . Thus an accurate determination of the amount of Pt 
in the C-layers is difficult from the reflectivity measure- 
ment. However, with the x-ray standing wave method 
it is possible to determine the amount of dissolved Pt 
in the C-layers through the detection of its fluorescence. 
Here the detection of Pt is direct and the fluorescence 
yield variation with angle for Pt in the C-layers has the 
opposite trend compared to Pt in the Pt-layers (see Fig. 
0) . So an analysis of the shape of the measured Pt fluo- 
rescence yield curve can provide the amount of dissolved 
Pt in C. 

We assume the presence of some dissolved Pt in C. 
Out of total Pt a fraction f c of Pt remains in the Pt- 
layers and within the broadened interface regions of the 
C-layers, and the remaining fraction (1 — f c ) is dissolved 
uniformly in the C-layers. Pt fluorescence yields as a 
function of angle for / c =l, 0.9 and 0.8 are shown in Fig. 
^. Later we will show with the experimental data that the 
fit to fluorescence yield improves when an / C <1 is allowed 
in the least-squares fitting procedure. From the amount 
(1 — f c ) of Pt in the C-layers we can obtain the average 
composition (Pt^ C±- x ) of the C-layers. In the present 
example, fractions / c =l, 0.9 and 0.8 correspond to 0, 4.4 
and 8.7% Pt in the C-layers, respectively. Keeping / c =l, 
it is also possible to fit the fluorescence data assuming 
broader interfaces, i.e., allowing larger values of o\ and 
(72 in Eqns.(28) and (29). However, it would be incon- 
sistent with the values of o~\ and oi obtained from the 
analysis of reflectivity data, as will be shown in section 
IV. In order to obtain a consistent set of microstructural 
parameters, it is necessary to allow, an f c < 1. f c may be 
called coherent fraction and (1 — f c ) incoherent fraction 
in analogy with the XSW analysis with Bragg diffraction 
from single crystals 41 1. 



III. EXPERIMENTAL 



Pt/C periodic multilayers with different bilayer pe- 
riod lengths d ranging from 35 A to 47 A were made on 
float glass substrates, kept at room temperature, by dc 
magnetron sputtering specially designed for coating in- 
ner walls of cyllindrical surfaces. Two sputter sources of 
Pt and C are located at top and bottom of the cyllindri- 
cal vacuum chamber. Samples were grown at low argon 
pressure of 1 mbar. The deposition rate of Pt and C was 
1 A/sec and 0.4 A/sec, respectively. The layer thick- 
ness during deposition was controlled using ion current 
and deposition time. Uniformity in the horizotal plane is 
achieved by rotating the sample while vertical uniformity 



is acieved by the mask. The overall thickness variation 
was found to be <2% over an area of 10 cmx 10 cm. The 
control of the thickness of individual layers was within 1 
A. A total of 20 layer pairs of Pt/C were deposited in 
each case. The x-ray specular reflectivity measurements 
have been made on these samples [Q to determine the 
bilayer thickness and interface roughness. We have used 
one of these samples for the combined x-ray standing 
wave and reflectometry analysis. 

Experiments were performed in our laboratory with a 
18 kW Mo rotating anode x-ray source. The experimen- 
tal set-up is schematically shown in Fig. |8[ Monochro- 
matic Mo K ai beam is obtained with the help of an 
asymmetrically cut Si(lll) crystal monochromator. The 
asymmetrically-cut crystal reduces the divergence of the 
monochromatized beam and is in standard use in X-ray 
standing wave experiments Q . The incident beam on the 
sample has an angular divergence of 0.006°. The vertical 
beam width is kept as small as 100 am. Reflected x-rays 
were detected with a Nal(Tl) detector and the Pt L a flu- 
orescent x-rays were detected with a Si(Li) detector. The 
reflected x-rays and the fluorescent x-rays were collected 
simultaneously at each angle. Control of the instruments 
for the operation of the HUBER diffractometer and data 
collection is obtained through a PC using Turbo C pro- 
gramming for IEEE and RS-232 protocols. More details 
about the set-up has been presented elsewhere J43| . The 
average exit angle a (the inclination of the Si(Li) detec- 
tor with respect to the sample surface) for fluorescent 
photons was 50°. 



IV. RESULTS AND DISCUSSIONS 

The experimental reflectivity data and the fitted theo- 
retical reflectivity curve (Theory-1) are shown in Fig. |^. 
Bragg peaks upto the third order are seen. The small 
oscillations are due to the total thickness of the multi- 
layer. Experimental data have been fitted by allowing 
the variation in the electron density, layer thickness, and 
surface and interface roughnesses of the layers. From the 
least-squares fitting the values of the parameters have 
been extracted. This fitting gives the Pt-layers density 
pi=4.95 electrons / A 3 , thickness d\— 16.8 A and C-layers 
density (fixed) p2=0.698 electron/ A 3 , thickness a\= 26.1 
A, A and 02=2.9 A. So the bilayer thickness is 

42.9 A. The third order peak position does not fit prop- 
erly. This may be due to the multilayer having a slight 
variation in bilayer thickness along the growth direction. 

It has been demonstrated that in case of single layer 
films the roughness is correlated with the thickness of 
the film But in the case of multilayer systems 

with alternating marker and spacer layers the roughness 
becomes complicated depending on the types of mate- 
rial, their diffusion properties, reaction and growth be- 
havior pC| |. It has been shown that in a W/C multilayer 
system the W-on-C interface is more rough than the C- 
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on-W interface p9| . Fundamentally this is expected be- 
cause of the nonwetting condition in the surface free en- 
ergy (aw > ad) f° r the growth of W on C. In our case 
apt > ac and we also observed the same trend: the 
Pt-on-C interface is more rough (<7i=4.5 A) than the C- 
on-Pt interface (<72=2.9 A). Pt electron density for this 
sample is 4.95 electron/ A 3 , which is lower than that of 
pure Pt electron density of 5.05 electron/ A 3 (p m =21.5 
gm/cc). In general, thin films tend to have a lower den- 
sity compared to pure bulk material. Additionally, in- 
tcrdiffusion across the interfaces leading to a mixed layer 
would decrease the Pt-layers density and increase the C- 
layers density. 

The Pt L a fluorescence yield has been measured over 
an angular region containing the first order Bragg peak 
and analyzed as follows. From the spectrum at each angle 
in the multichannel analyzer only Pt L a peak is selected. 
These peaks at all angles are fitted and the background- 
subtracted area is determined. This area gives the yield. 
This raw yield data have been corrected for footprint, 
probing thickness variation and finite detector aperture. 
These corrections are explained at the end of this section. 
This corrected Pt L a fluorescence yield vs. angle along 
with reflectivity over first Bragg peak is shown in Fig. 
|l0| . We fit the fluorescence yield data based on the model 
described earlier. This model incorporates all the param- 
eters extracted from the reflectivity fit. That means that 
the density, thickness, surface and interface roughness 
etc.. of the layers are kept intact. Here we have con- 
sidered the contribution of roughness as error functions 
[Eqns.(28) and (29)] at both the interfaces with cti=4.5 
A and (72= 2.9 A. These a- values are the roughness val- 
ues obtained from the analysis of reflectivity. (It is well 
known that reflectivity calculations using explicit error- 
function concentration profile at the interface and the 
flat interface reflection coefficient multiplied by a Debye- 
Waller function [Eqn. (7)] are equivalent [|45[ ). If we 
consider that there is no dissolved Pt in the C-layers (i.e. 
f c = 1), we do not obtain a good fit. The best fit is ob- 
tained with the model with a uniform mixing of Pt in 
the C-layers with / c =0.87. This means that 13% of total 
Pt is dissolved within the C-layers. Converted to atomic 
concentration this corresponds to the average composi- 
tion Pto.05Co.95 of the carbon layer. It should be noted 
that Pt concentration in the C-layers is actually higher 
near the interface. This concentration varies with the 
distance from the interface and can be easily determined 
from the distributions [Eqns.(28) and (29)]. 

In order to show the sensitivity of the fluorescence yield 
curve to the Pt concentration in the C-layers, we also 
show the plots for Pto.03Co.97 and Pto.07Co.93 in Fig. 
They are distinctly different from the data and the fit- 
ted curve for Pto.05Co.95. This clearly shows that the 
uncertainty in the estimated Pt concentration of 5% is 
smaller than 2%. In the fitting of data the weighted R- 
factors are 0.041, 0.031, 0.023, 0.024 and 0.029 for 3, 4, 
5, 6 and 7% Pt, respectively. It is noticed from Fig. 
that with increasing Pt concentration in C, Pt fluores- 



cence yield increases on the low- angle edge and decreases 
on the high-angle edge of the reflectivity curve. This can 
be easily understood from Fig. [| At angular position 
'b' on the reflectivity curve, X-ray intensity is high in 
the C-layers and low in the Pt-layers. However, if there 
is no Pt in the C-layers, there would be no Pt fluores- 
cence emission from there. As some Pt migrates from 
the Pt-layers to C-layers, the amount of Pt present in 
the C-layers would produce strong fluorescence emission. 
That is why increasing Pt concentration in the C-layers 
produces higher fluorescence yield at this angular posi- 
tion 'b' as seen in Fig. It is also noticed from Fig. || 
that the maximum field intensity in the C-layers is much 
higher than the maximum field intensity in the Pt-layers 
(see also Fig. |]). This is due to lower absorption of X- 
rays in the C-layers. Due to this fact, a given amount of 
Pt in the C-layers produces a stronger fluorescence signal 
than the same amount in the Pt-layers when the X-ray 
intensities are maximum in the respective layers. 

Probing the quantity of material dissolved from one 
layer into the other layer of a layer-pair in a multilayer 
system is not only important for optical mirrors and de- 
vices, but also very crucial for magnetic multilayers where 
interface broadening and alloying within the layers affect 
magnetic properties of multilayers. In magnetic multilay- 
ers with alternating layers of magnetic and non-magnetic 
materials, a small amount (even a few percent) of mag- 
netic impurity in the nonmagnetic layers can change the 
magnetic coupling and magnetoresistance. In fact, in 
magnetic multilayers with a wide range of Cui_ x Ni x 
(x=0.0A to 0.42) alloy spacer, the smallest amount of im- 
purity (:r= 0-04) has shown the largest change in magne- 
toresistance p9| . The magnetic impurity in the nonmag- 
netic layer of the multilayer may be an element other than 
the magnetic element present in the multilayer. Since x- 
ray fluorescence can identify the element, the distribution 
of such impurity elements in the multilayer can be deter- 
mined by XSW experiments |^8| . 

It must be mentioned here that the fluorescence data 
can also be fitted, without assuming the dissolved frac- 
tion (i.e. keeping / c =l), by allowing a\ and (72 to vary 
for the fluorescence fit. This fit is also shown in Fig. [lC| . 
However, the cr-values obtained from this fit ( 01 = 8.9 A, 
(72 = 4.2 A) are inconsistent with those obtained from 
the reflectivity fit. The computed reflectivity for these 
a- values, as shown in Fig. || (Theory-2), is very differ- 
ent from the measured reflectivity. This shows that this 
set of larger a- values does not represent correct inter- 
face roughness. This is probably the reason why a very 
large a- value (10 A) fitted the fluorescence data of Kawa- 
mura et al. |24|]. Our results underline the necessity for 
the combined x-ray standing wave and reflectivity anal- 
ysis of periodic multilayers. We suggest that a combined 
use of reflectivity and x-ray standing waves can provide 
the microstructural details of a periodic multilayer. The 
procedure to follow is: (i) obtain bilayer periodicity, frac- 
tional thickness of the high-Z layer and surface and in- 
terface roughnesses from the reflectivity fit, (ii) interface 
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roughness should not be constrained to be equal for both 
types of interfaces, and (iii) use the parameters obtained 
from the reflectivity fit and proceed for the the fluores- 
cence data fit with the assumption of a dissolved fraction 
of one material in the other, either in uniform distri- 
bution or with any other improved distribution model. 
For a more accurate determination of this distribution, 
higher order Fourier components of the distribution can 
be determined by XSW measurements with higher order 
Bragg diffractions. 

In order to fit the reflectivity data to Eqn.(13) and 
fluorescence data to Eqn.(27), the following corrections 
to data were applied: (i) Footprint correction was 
applied to both reflectivity and fluorescence data. At 
very small angles the beam projection is larger than the 
sample area. So, only a fraction of incident photons are 
actually incident on the sample. After this correction, 
the data represent what they should be if all the photons 
were incident on the sample, (ii) The fluorescence data 
come from a relatively thin layer (thickness of the mul- 
tilayer) compared to the beam penetration depth. Thus 
with the variation of 8 the effective probe depth changes. 
To correct for that, fluorescence data are to be multiplied 
by sind at each point, (iii) The fluorescence detector has 
a finite aperture and the fluorescent photons may come 
from a much larger sample area. The detector offers a 
varying effective solid angle for fluorescent photons orig- 
inating from different parts of the sample surface. As 
the exposed sample area varies with 8, this requires a 
correction which depends on detector aperture, detec- 
tor distance from the sample and the sample length. In 
our case, over the 8 range (0.45°— 0.6°) of the first order 
Bragg peak region this introduces only a minor correction 
for 1% variation in detected intensity. 



V. CONCLUSIONS 

For a periodic multilayer system with alternating lay- 
ers of a high-Z and a low-Z element, Bragg diffraction of 
x-rays occurs when the Bragg condition for the bilayer 
periodicity is satisfied. As in diffraction from a large 
perfect crystal, standing waves are set up in the multi- 
layer while diffraction occurs. The antinodal (or nodal) 
planes of the standing wave are parallel to the layer- 
planes and have the periodicity equal to the multilayer 
period. On the low-angle side of the Bragg-reflection 
peak, the antinodal planes are within the layers with low- 
Z element. As the angle of incidence advances through 
the diffraction peak the antinodal planes shift inward and 
finally coincide with the nearest layer of high-Z element 
of the layer-pairs. Emission processes, such as photoe- 
mission or fluorescence from atoms in the multilayer are 
modulated, over an angular region containing the Bragg 
peak, following the shift of the antinodal planes. Analysis 
of this modulation in the emission yield provides struc- 
tural information about the multilayer. The usefulness of 



the combined application of x-ray reflectivity and x-ray 
standing wave techniques for the analysis of multilayer 
microstructures has been explained. Deficiencies of each 
technique can be overcome by the combined application 
of these techniques. XRR depends on the electron den- 
sity difference between the layers of the bilayer. Where 
electron density of one layer of the layer-pair is very small 
compared to the other, reflectivity is not very sensitive 
to even a large fractional change of this electron density. 
Moreover the change of electron density is not necessar- 
ily due to the diffusion of atoms from the other layer of 
the layer-pair, it could also be due to other impurities 
incorporated during multilayer fabrication. Thus accu- 
rate determination of the layer composition from XRR 
technique is practically impossible. These aspects have 
been elucidated with an example of a 20 period Pt/C 
multilayer. In the XSW technique, elements are directly 
identified. Thus the amount of dissolved Pt or any other 
impurity in the C-layers, such as Ar, often incorporated 
during multilayer fabrication, can be determined. As 
interface roughness drastically affects the higher order 
Bragg peaks and overall intensity at higher angles, in- 
terface roughnesses are more accurately determined by 
fitting the reflectivity data over a large range of angle 
of incidence. On the other hand, in the XSW analy- 
sis, if the amount of Pt in the C-layers is assumed to 
be solely within the broadened interface and treated as 
roughness, one obtains too large roughness values com- 
pared to those obtained from the reflctivity fit. Fixing 
the interface roughness values at those obtained from the 
XRR analysis and assuming the remaining Pt to be in 
uniform distribution in the C-layers, the Pt concentration 
in the C-layers is determined. (More details about the 
elemental distribution, such as higher order Fourier com- 
ponents, can be obtained by XSW measurements with 
higher order Bragg peaks). Thus a combined analysis by 
XSW and XRR techniques removes the deficiencies of the 
individual techniques. For a 20 period Pt/C multilayer 
system interface roughnesses ( Pt-on-C: 4.5 A, C-on-Pt : 
2.9 A) and the C-layers composition ( Pto.05Co.95) have 
been determined. Determination of a small quantity of 
impurity, even a few percent, in the spacer layer is par- 
ticularly important in the magnetic multilayers. 
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FIG. 1. A schematic representation of x-ray reflection 
from a multilayer system. See text for details. 
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FIG. 2. Reflectivity from a 20 period Pt/C multilayer sys- 
tem with periodicity d(43A) = di(17A) + d 2 (26A) and with 
surface and interface roughnesses (A) oo, <Ti, 02 0, 0, 3, 3, 

3 ( ) and 3, 5, 3 0, 0, ( ) 3, 3, 3 ( ) and 3, 5, 

3 ( ). 
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FIG. 3. X-ray standing wave field intensity distribution 
within the Pt/C multilayer system, at different angles of in- 
cidence 9 over the 1st order Bragg peak region (shown in the 
inset), (a) 9 = 0.35° ( ), (b) 9 = 0.486° ( ), (c) 



-), (e) 9 = 0.535° (- 



and 



9 = 0.500° ( ), (d) 9 = 0.516° (- 

). The phases, v(9), of E-field ratios ^ ( 

( ) are also shown in the second inset, which also shows 

reflectivity ( , x3). At a given depth z, the variation 

in field intensity with angle over the strong reflection region 
occurs mainly because of large variation in phase, v(9). [ see 
Eqn.(22)]. 
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FIG. 4. Integrated XSW field intensity over the Pt-layers 
and over the C-layers for different surface and interface rough- 
nesses for the Pt/C multilayer system. Reflectivity over the 
1st order Bragg peak (solid squares) for ao—O, o"i=0, o"2=0 
(A), integrated field intensity over Pt-layers with cro,c"i,o"2 

(in A): 0, 0, ( ); 3, 3, 3 ( ); 3, 5, 3 ( ); 3, 5, 

5 ( — ); 3, 7, 7 (• • •) and integrated field intensity over 

C-layers (connected open circles) for o"o=3, (Ji=5, <72=3 (A). 
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FIG. 5. Schematic diagram of Pt distribution, f(z), with 
interface roughness over the bilayer period. Co, <J\ and 02 
are surface roughness, Pt-on-C and C-on-Pt interface rough- 
nesses, respectively 
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FIG. 7. Theoretical plots of reflectivity for different elec- 
tron densities pc of the C-layers. pc =0.698 electrons/ A 3 
( ) t ^0=0.803 electrons / A 3 (15% higher compared to 

the actual density) ( ). Curves are vertically shifted by 

two orders. However, they are also shown in overlaping mode 
to demonstrate that they are practically indistinguishable. 
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FIG. 6. Theoretical plots for Pt fluorescence yield, com- 
puted for the distribution of Pt in Fig. ^, over the first order 
Bragg reflection angular region. Reflectivity (solid squares), 
Pt fluorescence yield integrated over Pt-layers with surface 

and interface roughnesses A, ai=5 A, <T2=3 A ( ), 

Pt fluorescence yield integrated over the whole multilayer 

(<T = 3 A, (Ti= 5 A, 0-2= 3 A and for / c =l ( ), / e =0.9 

( ) and / c =0.8 ( ). See text for details. 




FIG. 8. A schematic view of the experimental set up with 
an asymmetric Si(lll) crystal monochromator (MC) and in- 
cident x-rays from an 18 kW rotating Mo anode x-ray gener- 
ator. Slits: So, Si (horizontal width = 4 mm, vertical width 
= 100 pm), S2 (horizontal width = 10 mm, vertical width 
= 150 pm); Dl : Nal(Tl) scintillation detector ; D2 : Si(Li) 
energy dispersive detector; S: sample. 
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FIG. 9. Experimental reflectivity data (o o o) and fitted 

theoretical reflectivity curve ( ) for a Pt/C multilayer 

on a glass substrate with 20 bilayers. Parameters obtained 
from the fit: bilayer thickness d = 42.9 A, Pt layer thick- 
ness di = 16.8 A and C-layers thickness t?2 = 26.1 A, surface 
roughness oo = 3 A, Pt-on-C interface roughness o\ = 4.5 A 
and C-on-Pt interface roughness 02 = 2.9 A. Theoretical 

reflectivity curve ( ) for oo = 3 A, o\ = 8.9 A and 

<T2 =4.2 A and all other parameters unchanged. See text for 
details. 
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FIG. 10. Experimental Pt L a fluorescence yield ( O O O 
) and reflectivity (solid squares) vs angle of incidence 8 over 
the first order Bragg reflection and the theoretical curves: 

( ) (To=3 A, <7i=4.5 A, <72= 2.9 A, / c =1.0 (no Pt in 

C-layers); ( ) o =3 A, oi=4.5 A, o 2 =2.9 A, / c =0.87 

(Pto.05Co.95); ( ) <To= 3 A, (Ti-8.9 A, (72=4.2 A, /c = 1.0. 

Also shown, (for o =3 A, ai= 4.5 A and 02=2.9 A) are the 
curves (open squares) for Pto.03Co.97 (/ c =0.935) and (filled 
circles) for Pto.07Co.93 (/ c =0.844). Fluorescence curves have 
been normalized at #=0.4° as in Fig. rj. See text for details. 



12 



